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SUMMARY 


A technique to measure remotely sea-surface temperature and salinity was 
demonstrated with a dual-frequency microwave radiometer system. Accuracies in 
temperature of C and in salinity of 1 part per thousand (®/oo) for salinity 
greater than 5 ®/oo were attained after correcting for the influence of extra- 
terrestrial background radiation, atmospheric radiation and attenuation, sea- 
surface roughness, and antenna beamwidth. The radiometers, operating at 1.43 
and 2.65 GHz, comprise a third-generation system using null balancing and feed- 
back noise injection. Flight measurements from an aircraft at an altitude of 
1.4 km over the lower Chesapeake Bay and coastal areas of the Atlantic Ocean 
resulted in contour maps of sea-surface temperature and salinity with a spatial 
resolution of 0.5 km. 


INTRODUCTION 

Measurements of ocean temperature and salinity are useful for studying the 
circulation in bay areas and tracing river outflow. In particular, estuaries, 
as described by Thomann (ref. 1), are areas of considerable dynamic change in 
water temperature and salinity, being areas of mixing between fresh and ocean 
water. A disturbance of the status quo, caused for example by hurricanes, 
droughts, or extreme hot or cold spells, can cause severe stress on estuarine 
inhabitants (ref. 1). Also, maps of ocean temperature are useful for tracing 
current systems such as the Gulf Stream. Some ocean fish tend to swim along 
oceanic isotherms at certain times, so that a knowledge of temperature can 
assist in their location. 

Because the microwave emission of ocean water is dependent on temperature 
and salinity (ref. 1), it is possible to utilize microwave radiometers to sense 
sea-surface temperature and salinity remotely. The emission is more dependent 
on salinity at lower microwave frequencies and more dependent on temperature at 
higher microwave frequencies. Thus, a system comprising two microwave radiome- 
ters operating at 1.43 and 2.65 GHz can provide remote measurement of both param- 
eters. The use of a dual-frequency microwave radiometer system has an advantage 
over other types of remote sensing instruments in that it can operate at night 
and probe through clouds. Previously reported (ref. 2) measurements of salinity, 
obtained with a single-frequency microwave radiometer, had a root-mean- square 
accuracy of 2 to 3 °/oo; however, low-altitude measurements with an infrared 
radiometer were required to obtain sea-surface temperature. 

In this report, a dual-frequency microwave radiometer system is described 
along with the necessary corrections to the radiometer apparent temperature for 
various atmospheric and sea-surface effects- Although these corrections are 
small, they are necessary for accurate inversion of the measurements to obtain 
sea-surface temperature and salinity to the desired accuracy of 1® C and 1 ®/oo. 



Inversion algorithms developed for this purpose were used to produce contour 
maps of these quantities. 

SYMBOLS 

e emissivity at nadir of calm sea surface, dimensionless 

e(6,Tg) emissivity of sea at molecular temperature Tg for viewing angle 0, 
dimensionless 

ej^ emissivity at wavelength X, dimensionless 

F values of sea-surface temperature or salinity 

f frequency, GHz 

h aircraft measurement altitude, km 

P((|>) antenna power pattern, dimensionless 

S salinity, parts per thousand (®/oo) 

<T> average physical temperature of intervening atmosphere, K 

Tg equivalent antenna temperature, K 

"^atm terrestrial sky brightness temperature at viewing angle 0, K 

Tg sea-surface brightness temperature, K 

Tq total extraterrestrial background radiation temperature, K 

Tcos thermal component of cosmic background radiation, K 

Tgai galactic noise temperature, K 

TgCh) equivalent atmospheric noise temperature at altitude h, K 
Tr radiometer apparent temperature, K 

(Tp) average apparent temperature, K 

Tg sea-surface molecular temperature, K 

T(z,0) atmospheric physical temperature at altitudes z and viewing 
angle 0 , K 

^"^cor total correction to radiometer apparent temperature, K 

ATp antenna pattern temperature correction, K 
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sea-surface roughness temperature correction, K 
w wind speed, m/sec 

regression coefficients (see eq. (23)) 
z altitude, km 

a atmospheric absorption coefficient, dimensionless 

h(<|>i) antenna beam efficiency at 4>-| , dimensionless 
0 viewing angle from nadir, rad 

X wavelength, cm 

T(h) atmospheric opacity at altitude h, dimensionless 

Tq total one-way atmospheric opacity, dimensionless 

t(z) atmospheric opacity at altitude z, dimensionless 

<ti antenna pattern angle from nadir, rad 

4)1 antenna half-beamwidth for 98-percent beam efficiency, rad 

Superscripts: 

L 1.43 GHz or L-band 

S 2.65 GHz or S-band 


DESCRIPTION OF RADIOMETER SYSTEM 

The L-band (1.43 GHz) and S-band (2.65 GHz) radiometers used in this inves- 
tigation are third-generation, advanced, switched radiometers. The design of 
these radiometers is based on the application of two concepts (ref. 3) to the 
switched-input-type receiver introduced by Dicke (ref. 4). The first concept 
is to equalize the temperature of the reference noise source at the circulator 
switch with the temperature of the lossy microwave components between the antenna 
terminal and the receiver input. This temperature is maintained extremely con- 
stant (±0.03 K) . The second concept is to use a feedback loop to inject pulsed 
portions of a constant noise source (avalanche diode) into the received noise 
power at the antenna terminals to equalize the noise power at both inputs of the 
circulator switch. The pulse frequency, which determines the average value of 
injected noise power, is a measure of the noise power (radiation) received by 
the antenna. As a result of these two design improvements, the radiometers are 
nearly independent of gain variations and errors that are contributed by front- 
end losses. The radiometers therefore exhibit the long-term stability that is 
necessary to achieve absolute brightness temperature measurements to within a 
few tenths kelvin. The stability also eliminates the need for periodic calibra- 


3 



tion cycles and the resulting loss of data. The S-band radiometer, described in 
detail in reference 3» has a superheterodyne receiver. The L-band radiometer, 
developed at the NASA Langley Research Center, is similar to the S-band radiome- 
ter except that the L-band radiometer has a direct-type receiver; that is, it 
has no mixer and intermediate-frequency (IF) section.. 

A simplified block diagram of the radiometer system is shown in figure 1 . 

The radiation received by the antenna has the injected noise added to it. The 
circulator switch then switches between this sum of noise power and that of the 
matched termination which is at the reference temperature. The circulator 
switch output is amplified by the receiver's low-noise parametric amplifier and 
detected by a square-law detector. The detector output is demodulated in syn- 
chronism with the circulator switch to produce a voltage driving the voltage- 
to-frequency converter. The output signal from the voltage- to- frequency conver- 
ter is used to gate the noise injection pulse driver which adjusts the injected 
noise power level to achieve a null at the output. Gating the noise source with 
constant-width pulses produces a linear relationship between the pulse frequency 
and the average injected noise. It is this linearity that allows accurate mea- 
surement over a wide range of brightness temperatures (typically 303 K). Finally 
the output data are recorded in digital form on magnetic tape and in analogue 
form on a strip-chart recorder. In both cases, even slight variations in the 
reference temperature are recorded for correction of the output data. The 
S-band and L-band radiometers have resolutions of ±0.08 K and ±0.09 K, respec- 
tively, and with suitable calibration, both have an absolute accuracy of better 
than ±0.2 K. 

The output data of both radiometers are converted to digital form by a data 
processor developed at the NASA Langley Research Center. The processor also 
conditions and formats the housekeeping data from other sources that are neces- 
sary for the reduction of the radiometric data, such as flight parameters, time, 
latitude, and longitude. The data recorded on digital magnetic tape are also 
displayed as L-band and S-band radiometer apparent temperature by the data pro- 
cessor. The data processor is capable of adjusting measurement integration times 
independent of the radiometer settings. This capability provides an efficient 
way to adapt the overall integration time to the aircraft altitude and measure- 
ment spatial resolution' (antenna half-power footprint size). 


MICROWAVE BRIGHTNESS TEMPERATURE OF THE SEA 

Quantitative measurements of the thermal emission from the sea surface 
require corrections for the radiative and transmissive properties of the inter- 
vening atmosphere and sky, antenna beam efficiency, and sea-surface roughness. 
At the measurement wavelengths of 11 and 21 cm, there are four principal atmo- 
spheric and sky effects which must be considered: 

(1) Cosmic background radiation 

(2) Radiation from galactic and discrete celestial radio sources 
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(3) Attenuation from oxygen and water vapor 

(4) Radiation from the Sun 

It has long been known that the Earth's atmosphere is essentially transpar- 
ent to transmission of electromagnetic radiation at frequencies of 1 to 3 GHz. 
Extensive work over the years on microwave signal propagation through the atmo- 
sphere at centimeter wavelengths has indicated that the influence of clouds is 
small at these frequencies except under very severe storm conditions. An added 
factor for consideration is that the background galactic noise tends to decrease 
substantially as frequencies increase beyond about 1 GHz. Therefore, the fre- 
quency regime from 1 to 3 GHz is a well-suited choice for minimizing the effects 
of extraterrestrial background radiation and atmospheric interference. 

Despite these advantages, accurate surface temperature measurement by air- 
borne radiometers in this microwave region requires detailed knowledge of these 
effects for correcting the instrumental observations. The corrections to the 
measured brightness temperature of the ocean surface can still be of the order 
of a few kelvins and therefore, must be taken into account. The apparent tem- 
perature Tp (which may also be called the equivalent radiometric temperature 
of the complete set of received radiations) is calculated from the equation of 
radiative transfer by making use of the Rayleigh-Jeans approximation to the 
Planck law: 


Tr(0) = {l - e(0,Ts)).<|Tc exp [-/o“ a(0,z) dzj 

T(z,0) exp a(0,z) dzja(0,z) dz^ exp a(0,z) dzj 


f 


+ e(0,Ts)Tg exp 


-J^ a(0,z) dz + I T(z,0) exp 

J Jq 


1^-^^ a( 0 ,z) dz 


a(0,z) dz j a(0,z) dz 

( 1 ) 


where 


e(®>Tg) emissivity of sea at molecular temperature Tg for nadir viewing 
angle 0 

h altitude of aircraft 

Tq total extraterrestrial background radiation temperature 

T(z,0) temperature of atmosphere at altitudes z and viewing angle 0 
a(0,z) atmospheric absorption coefficient at z and 0 

The first term in equation (1) comprises the temperature of the downward radia- 
tion of the extraterrestrial noise, attenuated by the entire atmosphere, and the 
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downward radiation of the atmosphere itself, reflected by the ocean surface and 
in turn attenuated by the intervening atmosphere between the ocean and radiome- 
ter. The second term accounts for the attenuated emission from the ocean sur- 
face. ■ The third term accounts for the upward emission of the atmosphere between 
the ocean and the radiometer. 

Since the antenna gain pattern for most radiometers is highly directional 
and falls off rapidly with increasing 6 , a horizontally stratified atmosphere 
model can be used to approximate equation ( 1 ) by replacing the nested integrals 
with opacities: 

Tr( 0) = (l - e(e,Tg)}|TQ exp [-Tq sec 0^ 

+ T(z) exp [-t(z) sec 0]a(z) sec 0 dzj> exp [[-T(h) sec 0^ 

+ e(0,Tg)Tg exp |^-T(h) sec 0^ + T(z) exp [^-t(z) sec 0^a(z) sec 0 dz 

( 2 ) 

where Tq is the total one-way opacity of the atmosphere and T(h) and t(z) 
are the opacities at altitudes h and z. 

The extraterrestrial background radiation (Tq) consists of galactic noise 
(Tgai)> the thermal component of cosmic background radiation (TgQg), and emis- 
sion from strong radio sources such as the Sun. Since the Sun is a point source, 
its contribution depends strongly on the antenna viewing angle, the Sun's posi- 
tion, the antenna pattern, and the roughness of the sea surface. However, as 
explained subsequently, the Sun's contribution was considered to be negligible. 
Therefore, 


To = Tgog + Tga2 (3) 

The galactic radiation component includes only that portion of the galactic 
radiation entering vertically into the measurement area. The galactic sector 
6 hours to 9 hours right ascension and 27° to 47° declination was overhead dur- 
ing the time of measurement and within the 20° antenna beam coverage. The ther- 
mal contribution of the galaxy is determined by a curve-fitting method suggested 
in reference 5. The equivalent exo-atmospheric temperature in kelvins of the 
galactic radiation from the specified galactic sector during the measurement 
period is 


Tgai = 2.34f-2*53 (4) 

where f is the frequency in GHz. The cosmic background radiation TgQg is 
frequency independent over the microwave band (ref. 6) and is equal to 2.7 K. 

For the clear atmosphere, water vapor and molecular oxygen are the major 
sources of attenuation in the microwave region. At frequencies of 1.43 and 
2.65 GHz, attenuation due to water vapor is about an order of magnitude less 
than that due to molecular oxygen. Since the total opacity of the atmosphere 
is only of the order of 0.01, it is sufficient to consider the absorption due 
to molecular oxygen alone. The brightness temperature Tgtm the terrestrial 
sky at zenith angle 0 is 
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(5) 


T(z) exp |]-t(z) sec eJaCz) sec 0 dz 

which contains the atmospheric term of equation (2). It has been measured for 
the clear atmosphere at 2.65 GHz (ref. 7). These observations indicated that 
the zenith brightness temperature and, therefore, the opacity of the cloudless 
terrestrial sky, are consistent with the predictions of the Van Vleck-Weisskopf 
model for oxygen absorption. The measurements at 2.65 GHz give 2.2 K for the 
zenith sky temperature and a value for Tq of 0.0091. For 1.43 GHz, the theory 
predicts values of 2.1 K and 0.008 for T^tm "^o> respectively. The term 

e{0,Tg)Tg of equation (2) may be replaced by the sea-surface brightness temper- 
ature Tg. 

Because Tq and x(h) are very small, the exponential terms of equa- 
tion (2) may be approximated by 

exp [^-T sec 0^ = 1 - T sec 0 (6) 

Also, by considering that the measurements of this investigation were taken 
along the nadir (i.e., 0 = 0), sec 0 in equation (6) can be set to unity. 

A further simplification is made by evaluating the one-way opacity between the 
sea surface and the radiometer antenna as 

T(h) = a(z) dz (7) 

and the equivalent noise temperature of the intervening atmosphere as 

To(h) = T(z) exp [-T(z)]a(z) dz (8) 

Equation (8) is the last term of equation (2) with sec 0=1. The atmospheric 
opacity or attenuation T(h) and equivalent noise temperature TgCh) of the 
1962 U.S. Standard Atmosphere at 1.43 and 2.65 GHz for different altitudes are 
shown in figure 2. Figure 2 also shows that at altitudes less than 2 km, the 
opacity x(h) is approximately linear with altitude, having values of approxi- 
mately 0.00136h and 0.00154h at 1.43 and 2.65 GHz, respectively. The right 
side of equation (8) can (for h < 2.5 km) be approximated by <T> T(h) where 
<T> is the averaged physical temperature of the intervening atmosphere between 
the radiometer and the sea surface. 

By applying the preceding approximations and derivations and replacing the 
term e(0,Tg) with e (the nadir emissivity of the calm sea surface), equa- 
tion (2) simplifies to 

Tr = TbD - T(h)] + (1 - e)[l - T(h)] RTcos + ‘ + Tatml + ”r(h) <T> 

(9) 

Unless the sea surface is perfectly smooth (no waves), a surface roughness cor- 
rection (AT^,) is necessary at some microwave frequencies because of the depen- 
dence of emissivity on roughness. For the exact measurement of the apparent 
temperature, the antenna should have a "pencil” beam pattern. Since this is 
impossible to achieve, an antenna pattern correction (ATp) must be made (ref. 8). 
Adding these two additional corrections to equation (9) results in 
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( 10 ) 


Tr = Tb[i - T(h)] + (1 - e) [l - T(h)][(Tcos + “ '^o> + 

+ T(h) <T> + At„ + ATp 

The sea-surface roughness can be expressed in terms of wind speed w. On 
the basis of prior flight data obtained with the S-band radiometer (ref. 9), 
the wind speed correction at 2.65 GHz was found to be 0.56w®*53. n© correction 
for surface roughness due to wind speed at 1.43 GHz is included because recent 
data indicate that it is negligible at this frequency (ref. 10). 

The antenna pattern correction is obtained by evaluation of the antenna 
temperature integrals (ref. 9), 

Tr( 4>) P((|>) sin 4> d4> TR((j>) P(4)) sin 4> d(j> 

P(4>) sin 4> d(l) P(4)) sin (j) d<l> 

where the antenna temperature T^ is the equivalent temperature of the radia- 
tion incident within the antenna beamwidth (24>i) plus the radiation incident 
outside this beamwidth, Tr(<|>) is the apparent temperature of the sea surface 
in the direction ()> from the nadir, and P(<j>) is the antenna power pattern. 

The antenna beamwidth 2<t>i for a 98-percent antenna beam efficiency was 50° in 
this case. The first term in equation (11) is by definition the antenna beam 
efficiency T)(4 >i). For antennas with high beam efficiencies, such as those of 
the L- and S-band radiometer system, equation (11) may be approximated by 

Ta = 0(<l>i) Tr + [l - ti(<|)i)]<Tr> (12) 

where (Tr) is the average apparent temperature from radiation present outside 
the antenna beamwidth. The antenna pattern correction is then 

^Tp = Ta - Tr = [l - T1 (<{>i)][<Tr> - Tr] (13) 

Evaluation of equation (13) for the 2.65-GHz antenna is described in detail in 
reference 9 and was found to yield 0.4 K. A similar evaluation was performed 
with the aid of the 1.43-GHz antenna patterns shown in reference 11 and resulted 
in an antenna pattern correction of 0.14 K. 

The results of evaluation of all unwanted contributions to the observed 
temperature (radiometer apparent temperature) are listed in table I for both 
frequencies. Using these results in equation (10) gives 

TrS = Tb^[i - xS(h)] + [l - eS][l - xS(h)]5.1 + TS(h) (T> + 0.56w0-53 + 0.4 

(14) 

TrL = Tb'^[i - TL(h)]| + [l - e^][l - xL(h)]5.8 + xL(h) (T> + 0.14 (15) 

where the superscripts S and L refer to 2.65 GHz and 1.43 GHz, respectively, 
TrS and Tr^ are the respective radiometer apparent temperatures, h is the 
aircraft altitude in kilometers, and w the wind speed in meters/second. Equa- 
tions (14) and (15) may be used without affecting the measurement accuracy if 
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(1) measurement is in nadir direction, (2) there is no precipitation, (3) the 
radiometers are operated at microwave frequencies lower than 3 GHz, (4) the 
altitude of the radiometers is low (less than 2.5 km), and (5) the Sun is low 
in the sky (i.e., negligible Sun contribution). 

Rearranging equations (14) and (15) gives 

TrS = [l - xS(h)] [tb^ + 5.1(1 - eS)J + xS(h) <T> + 0.56w0-53 + 0.4 (16) 

TrL = [l - TL(h)][TB^ + 5.8(1 - eL)] + xL(h) <T> +0.14 (17) 

Because the terms 5.1(1 - e^) and 5.8(1 - e^) are always small in comparison 
with Tb^ and Tb^, the (1 - e) terms may be approximated by the constant 
value of 0.65 since the emissivity e is close to 0.35. Hence, 

= Tb^ + 3.7 + xS(h)[<T> - Tb^ - 3 . 3 ] + 0.56w0-53 (18) 

TrL = Tb^ + 3.9 + xl^(h)[<T> - Tb^ - 3.8] (19) 

Within the bracket term only, <T> , Tb^, and Tb^ may normally be approximated 
by 283 K, 105 K, and 95 K, respectively. 

Using the expressions listed in table I for xS(h) and x^(h) and solving 
for the brightness temperature of the sea surface Tb result in the following 
two equations; 


TbS = 

TrS _ 3.7 _ o.269h - 0.56w°-53 

( 20 ) 

Tfi^ = 

TrI' - 3.9 - 0.251h 

( 21 ) 


DETERMINATION OF SALINITY AND SURFACE TEMPERATURE 
FROM MICROWAVE BRIGHTNESS TEMPERATURE 

The brightness temperature Tb measured by an ideal radiometer is related 
to the molecular temperature of a radiating surface via the emissivity of the 
surface. The emissivity of a dielectric surface at a particular wavelength is 
determined by its complex dielectric constant which for sea water is a function 
only of temperature and salinity. Therefore, the brightness temperature of the 
sea surface is given by 

Tb(X) = ex(Ts,S)Tg (22) 

where the emissivity e at the wavelength X is expressed in terms of surface 
temperature Tg and salinity S. Plots of brightness temperature as a function 
of salinity and surface temperature at 1.43 GHz and 2.65 GHz are given in fig- 
ures 3 and 4, respectively (ref. 12). These plots were basically derived from 
experimental measurements of the dielectric constant of sea water at 1.43 GHz 
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(ref. 13) with a few experimental errors removed by utilizing data obtained 
previously at 2.65 GHz (ref. 14). The frequency dependency was accounted for 
by means of the Debye equation (ref. 15). The values on which the family of 
curves are based have published accuracies (ref. 12) of 1 °/oo for salinity 
(S greater than 5 ®/oo) arid 0.5° C for sea-surface temperature. 

In order to invert the brightness temperatures Tg® and Tg^ in figures 3 
and 4 to corresponding values of Tg and S, the following regression equation 
was developed: 

F = XiTgS + X2Tg^^ + X3TgSTg^ + Xn(TgS)^ + X5(TgL)^ + X6(TgS)3 

+ X7(TgS)2TgL + X8 TbS(TbL)2 + X9(Tb^)^ (23) 

where F is equal to either Tg or S depending on the set of regression 
coefficients Xj^ (listed in table II) used. 

The regression coefficients were found from input data consisting of two 
9 by 9 tables of brightness temperature as a function of salinity and molecular 
surface temperature; one table was for 2.65 GHz and the other for 1.43 GHz. 

The coefficients X^ were computed so that the mean-square error between input 
salinity or temperature and the corresponding quantity computed from equa- 
tion (23) was minimized over the input data set. 


FLIGHT OBSERVATIONS AND COMPARISON WITH SEA-TRUTH DATA 

On August 24, 1976, the L- and S-band radiometer system was installed on 
a NASA C-54 aircraft, and a functional check-out flight was conducted from the 
NASA Wallops Flight Center over the lower part of the Chesapeake Bay and adja- 
cent Atlantic Ocean. A flight over the entrance to the Chesapeake Bay was 
selected because the mixing of fresh and salt water there should result in high 
salinity gradients, and instrument performance could therefore be verified over 
a wide range of salinity concentration. .The results are reported in the form 
of synoptic maps of salinity and temperature, and the data are compared with 
sea- truth measurements. These sea- truth measurements were obtained via the 
bucket method: water temperature of bucket samples taken at the surface was 

measured with a thermometer, and samples of the same water were analyzed in a 
laboratory for salinity using the silver nitrate technique. The accuracy of 
the sea-truth measurements was about 0.1 °/oo for salinity and about 0.1° C 
for water temperature. 

Sea conditions for the measurement were fairly calm with a 3.5-m/sec sur- 
face wind. During the test flight, the data were taken under cloudless but 
hazy atmospheric conditions with the aircraft in level flight (roll and pitch 
less than 2°) at an altitude of 1.4 km, so that all the measurements refer to 
observations at nadir. During the measurement time period between 9:00 a.m. 
and 12:00 p.m. EST, the elevation angle of the Sun varied from 41° to 64°, so 
that there was no Sun contribution to the nadir microwave measurements. The 
spatial resolution of the measurement, or antenna half-power footprint size, 
was 0.5 km. Introducing these parameters into equations (20) and (21) results 
in 
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5.2 


and 


where 


TfiS = TrS _ 

Tb^ a TrL - 
TgS and TgL. 


4.2 

are in ke Ivins. 


(24) 


(25) 


Both radiometers were calibrated before the flight by using a specially 
designed liquid-nitrogen cryogenic antenna load (ref. 16) placed under and 
closely coupled to the antenna. This calibration procedure is reported for 
the S-band radiometer in reference 9. Since all the corrections that were 
applied in deriving equations (24) and (25) were to some degree inaccurate 
and there were some inaccuracies associated with the radiometers in the mea- 
surement of microwave brightness temperature, an error budget has been prepared. 
Table III lists the possible errors in terms of radiometric brightness tempera- 
ture for each radiometer due to atmospheric and sky emission effects, antenna 
beam efficiency correction, calibration, and instrument stability. Included 
also is the temperature resolution (i.e., relative accuracy) for a 1 -second 
integration time and 100-MHz bandwidth for the S-band radiometer and 60-MHz 
bandwidth for the L-band radiometer. 


The measured radiometric data, together with latitude and longitude coordi- 
nates from the aircraft’s inertial navigation system, were recorded on digital 
magnetic tape and computer processed according to equations (23), (24), and (25). 
The sea-truth data were obtained within 30 minutes of the radiometric data from 
several locations in the measurement area. These locations and the flight lines 
used to obtain the radiometric data are shown in figure 5. Although it had been 
planned to have the sea-truth locations correspond to the flight lines, an accum- 
ulative drift of about 0.2 minute in latitude per flight line in the inertial 
navigation system during the measurement prevented this. A sample section of 
the raw analogue data taken along the first flight line (37^16’ latitude) and 
over sea-truth location 1 is shown in figure 6. The fluctuations of Tp^ and 
Tr^ around the mean values represent the fine structure of salinity and tempera- 
ture variations, since the brightness temperature resolutions of both radiometers 
are less than 0.1 K. It is obvious from the data traces that Tpj^ is much more 
responsive to salinity gradients between the Chesapeake Bay and the Atlantic 
Ocean than Tp^. The necessary corrections AT^or ^"^cor listed in table I 

and incorporated in equations (24) and (25) are indicated at sea-truth loca- 
tion 1 . The resulting brightness temperatures Tg^ and Tg^ are converted 
into a salinity of 18 °/oo and a temperature of 24.6^ C via the regression equa- 
tion (23). These values differ from the sea-truth surface measurement by less 
thain the target accuracy of 1 °/oo and 1° C. 

A comparison of the radiqmetrically measured values of surface temperature 
and salinity with those taken at the sea-truth locations is shown in figure 7. 
This figure shows a section of computer-derived remote measurements as a func- 
tion of latitude and longitude in the area of the surface measurements along 
with the values of the corresponding sea-truth data. Data points which were 
unusable because of radio interference from local radar or the presence of land 
are indicated by the symbol X. Several surface measurements were taken around 
sea-truth location 10 to determine variation of the sea-truth data in a given 
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small area. Since the flight lines did not always coincide with the sea- truth 
locations, a comparative analysis of the measured results with the sea-truth 
data required some arithmetic averaging of the measured data in the vicinity 
of the sea- truth locations. The values for salinity and surface temperature 
are compared with the corresponding sea- truth data in table IV. It should be 
noted that the sea- truth data shown for location 10 are mean values of the 
several measurements . Also , no comparison was made for location 11, as the radi- 
caneter measurements there were influenced by the presence of land. The result- 
ing error mean deviation of -0.52 ®/©o ahd 0.55° C and error standard deviation 
of 0.92 ®/po and 0.59° C for salinity and surface temperature, respectively, 
seem to indicate that the desired accuracy of 1 ®/oo salinity and 1° C for 
temperature was achieved, although only a 50-percent confidence factor can be 
applied to these deviations because of the small number of data points for 
comparison. 

The calculated values of salinity and surface temperature were plotted as 
a function of geographic position. A contour map of the salinity distribution 
is shown by isohalines in 2 ®/o© increments in figure 8. It can be seen from 
the figure that the higher salinities indeed exist in the ocean and are reduced 
considerably by mixing with fresh water outflow in the bay region. Evidently 
the ocean water with its higher salinities is pushed into the northern part of 
the Chesapeake Bay entrance by tidal effects and up the inside of the Delmarva 
Peninsula. The lower salinity bay waters flow south out the southern part of 
the bay entrance into the ocean, because of Coriolis forces (private communica- 
tion with B. Boicourt, CBI, Johns Hopkins Univ. , Baltimore, Md. , 1976). A con- 
tour map of the surface temperature in 2° C increments is shown in figure 9. 

The cooler temperatures exist in the ocean and are influenced by the warmer tem- 
peratures of the fresh water outflow in the region of mixing. The fine struc- 
ture of surface temperature caused by upwelling, heating in shallow waters by 
the Sun’s radiation, and outflow of rivers is averaged in order to show the 
general temperature distribution. Again, a southern outflow of the warmer 
Chesapeake Bay waters is indicated. 


CONCLUDING REMARKS 

A dual-frequency (1.^3 and 2.65 GHz) radiometer system was developed to 
measure salinity to 1 part per thousand (®/oo) (for salinity greater than 
5 ®/oo) and sea-surface temperature to 1° C. This system utilized a null- 
balancing feedback circuit with noise injection for stabilization. It was 
operated from an aircraft over the lower Chesapeake Bay and adjacent coastal 
area of the Atlantic Ocean to map salinity and sea-surface temperature. To 
obtain the stated accuracies, the radiometer data were corrected for extrater- 
restrial background radiation, atmospheric effects, sea-surface roughness, and 
antenna beam efficiency. 


Leingley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 26, 1977 
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TABLE I.- CORRECTIONS TO RADIOMETER APPARENT TEMPERATURE AT 1.43 AND 2.65 GHz 



Correction 

Comments 

X!il t 6C U 

2.65 GHz 
( S-band ) 

1.43 GHz 
(L-band) 

Cosmic background radiation, 
*^cos> ^ 

2.7 

2.7 

Summarized in reference 6 

Galactic radiation, 

K 

0.2 

1 .0 

Only for specified galactic 
quadrant according to 
equation (4) 

Atmospheric contribution 

(oxygen), Ta^m* K . . . . 

2.2 

2.1 

Consistent with Van Vleck- 
Weisskopf model and 1962 
U . S . Standard Atmosphere 

Nonideal antenna, ATp, 

K 

0.4 

0.14 

Both antennas are nonscan- 
ning , high-beam-efficiency , 
multi-mode circularly 
polarized horns with 20^ 
half-power beamwidth 

Surface roughness, AT„, 

K 

1.1 

0 

S-band correction obtained 
experimentally: 

AT„ a 0.56w°-5d 
L-band correction negligible 

Total atmospheric opacity, 

T^O 

0.0091 

0.008 


Opacity at altitude h, 

T(h) 

0.00154h 

0.00136h 

x(h) is linear up to 2.5 km 
(fig. 2) 

Total correction to Tp, 

^*^cor» ^ 

5.2 

4.2 

Cumulative correction factor 
to the apparent tempera- 
ture for w = 3-5 m/sec 
and h = 1 . 4 km 
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TABLE II.- REGRESSION COEFFICIENTS FOR SURFACE TEMPERATURE 


AND SALINITY COMPUTATION 



Regression coefficients for computation of 


Temperature 

Salinity 

Xl 

16.9073567947 

138.2129737430 

X2 

-21.8805724219 

-137.4748877279 

X3 

.4925788939 

7.0376869542 

X4 

-.3647221634 

-4.6052164921 

X 5 

-.0475842615 

-2.4460477714 

X6 

.0051508507 

.0403065628 

X 7 

-.0122197794 

-.0844284348 

X8 

.0099636042 

.0566040411 

X 9 

-.0031929574 

-.0124179422 


TABLE III.- ERROR BUDGET IN MEASUREMENT OF SEA-SURFACE 
BRIGHTNESS TEMPERATURE 


Error source 

S-band error, K 

L-band error 

Atmospheric and sky emission effects 

0.3 

0.3 

Antenna pattern 

.1 

.1 

Instrument calibration 

.1 

.1 

Instrument stability 

. 1 

.1 

Temperature resolution 

.08 

.09 

Maximum error, K 

0.68 

0.69 

Root-mean-square error, K 

.34 

.34 
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TABLE IV.- COMPARISON OF RADIOMETRIC AND SEA-TRUTH MEASUREMENTS 


Location no. 

Salinity, ^ /o 

O 

Temperature, 

Sea truth 

Radiometric 

Error^ 

Sea truth 

Radiometric 

Error^ 

1 

17.7 

18.0 

-0.3 

25.5 

24.6 

0.9 

2 

18.0 

20.0 

-2.2 

25.7 

25.5 

.2 

3 

19.0 

18.9 

.1 

25.7 

24.4 

+1.3 

4 

19.3 

19.9 

-.6 

26.0 

25.6 

.4 

5 

19.0 

20.4 

-1 .4 

25.6 

26.3 

-.7 

6 

18.4 

19.4 

-1 .0 

25.8 

25.9 

-.1 

7 

18.8 

18.5 

.3 

25.6 

24.3 

+ 1 .3 

8 

19.8 

19.6 

.2 

25.8 

25.2 

.6 

9 

20.8 

21 .8 

-1 .0 

26.3 

25.6 

.7 

10 

^21 .5 

22.4 

-.9 

t>26.2 

25.5 

.7 

12 

28.1 

27.0 

+ 1 • 1 

25.0 

24.3 

.7 

Error mean deviation 

Error standard deviation 

-0.52 °/oo 
.92 °/oo 


0.55° C 
.59° C 


®Error equals sea- truth value minus measured value. 
^Mean of value of several surface measurements . 
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Figure 2.- Atmospheric attenuation and equivalent noise temperature of the 1962 U.S. Standard 
Atmosphere for different altitudes at 1.43 GHz and 2.65 GHz. 
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Figure 5.- Flight pattern over the entrance to the Chesapeake Bay and sea-truth locations 
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Figure 7.- Comparison of remote measurements of salinity and temperature with surface measurements. 
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